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ABSTRACT 

Context. The HCN, HCO + , and HNC molecules are commonly used as tracers of dense star-forming gas in external galaxies, but 
such observations are spatially unresolved. Reliably inferring the properties of galactic nuclei and disks requires detailed studies of 
sources whose structure is spatially resolved. 

Aims. To understand the origin of extragalactic molecular line emission, we compare the spatial distributions and abundance ratios of 
HCN, HCO + , and HNC in W49A, the most massive and luminous star-forming region in the Galactic disk. 

Methods. Maps of a 2' (6.6 pc) field at 14" (0.83 pc) resolution of the J=4-3 transitions of HCN, H 13 CN, HC 15 N, HCO+, H 13 CO + , 
HC ls O + , and HNC are combined with supplementary observations of the J=5-4 transition of DCN and the J=3-2 transitions of HCN 
and H 13 CO + . Most of the data are from HARP/ACSIS, with supplementary data from JCMT Receiver A and the SCUBA archive. We 
use maps of the integrated intensity and line-profiles to pick out regions of the source to study in more detail. We compare column 
densities and abundance ratios towards these regions with each other and with predictions from gas-phase chemical models. 
Results. The kinematics of the molecular gas in W49A appears complex, with a mixture of infall and outflow motions. Both the 
line profiles and comparison of the main and rarer species show that the main species are optically thick. Two "clumps" of infalling 
gas that we look at in more detail appear to be at ~40 K, compared to > 100 K at the source centre, and may be ~10x denser than 
the rest of the outer cloud. The chemical modelling suggests that the HCN/HNC ratio probes the current gas temperature, while the 
HCN/HCO + ratio and the deuterium fractionation were set during an earlier, colder phase of evolution. 

Conclusions. The similarity in the derived physical conditions in W49A and those inferred for the molecular gas in external galaxies 
suggest W49A is an appropriate analogue of an extragalactic star forming region. Our data show that the use of HCN/HNC/HCO + 
line ratios as proxies for the abundance ratios is incorrect for W49A , suggesting that using these line ratios as abundance ratios in 
galactic nuclei is invalid too. On the other hand, our observed isotopic line ratios such as H 13 CN/H 13 CO + approach our modeled 
abundance ratios quite well in W49A. Second, the 4-3 lines of HCN and HCO + are much better tracers of the dense star-forming 
gas in W49A than the 1-0 lines, confirming similar indications for galactic nuclei. Finally, our observed HCN/HNC and HCN/HCO + 
ratios in W49A are inconsistent with homogeneous PDR or XDR models, indicating that irradiation does not strongly affect the gas 
chemistry in W49A. Overall, the W49A region appears to be a useful template for starburst galaxies. 

Key words. ISM: molecules - ISM: clouds - ISM: abundances - stars: formation 

1. Introduction trace the star forming activity as well as or even better than 

HCN 1-0 (iKru mholz & Thompso r32007tlNaravanan et alJ2008h 



H \ In recent years, submillimeter line emission from a variety 
of molecules has been det ected from external galaxies in the 



Bayet et al. 2009b). 



early and local Universe dKramer et al.l 120081) . The first de- The HCN / HCO + ratio is also used to trace the chemical 
tection of HCO + at high redshift, the J=l-0 transition to - state of the dense gas, in particular its irradiation by U V light 
wards the Cloverleaf quasar at z=2.56 (iRiechers et al.l 120061) . from a starburst or by X-rays from an AGN. For example, Kohno 
was quickly followed up by higher-/ HCN, HNC and CN (l2005h used the HCN/HCO + ratio in near by Seyfert galaxies 
lines tow ards the ultra-lum i nous quasar APM 08279+5255 at to distinguish between AGNs and starbursts : l(jreve et all J2009) 
z=3.91 1 ( Wag g et aUl2005t iGuelin et al.l 12007). Whereas low- find that excited HCN transitions trace dense gas in the two star- 
J CO lines are used to trace the total molecular gas mass in burst galaxies, Arp 220 and NGC 6240, while sub-thermally ex- 
Galactic nuclei, HCN and HCO + lines are often used to trace cited HCO + lines trace gas which i s 10-100 times le ss dense; 
the higher-density gas. This dense gas is much more closely re- the HCN/HCO + ratios observed by lAalto et alj d2007l) towards 
lated to star formation activity than CO, as sho wn by the tighter NGC 4418 suggest emission from dense, warm gas with a PDR 
corre lation of L(HCN 1-0) with L(far-IR) dGao & Solomon! component, but with no X-ray chemistry. Another commonly 
l200l . However, factors other than the mass of dense gas pre- used line ratio is HCN / HNC, whi ch is thought to trace the gas 
sumably affect the HCN emission, and high-./ CO lines may temperature (e.g jBaan et al. 2008). 
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The angular resolution of current instrumentation at submm 
wavelengths ranges from ~30" for single-dish telescopes to 
~\" for interferometers. For nearby galaxies (D ~10-100Mpc), 
single-dish observations thus probe scales of 1.5-15 kpc, while 
interferometers probe scales of 50-500 pc. In the case of high- 
redshift galaxies, even interferometers probe size scales of ~kpc. 
The goal of this work is to improve understanding of such ob- 
servations by using a local analogue. In particular, the use of 
line ratios to infer physical conditions depends on the assump- 
tion that the gas is roughly homogeneous on these scales, so that 
the lines trace the same gas. Recent authors have already ques- 
tioned this assumption (Bavetet aT] |2008l l2009al) . and the aim 
of this paper is to see just how far this assumption of homo- 
geneity is from reality, using the W49A region as a template for 
extragalactic regions with high star formation rates, and to ex- 
plore how the observational limitations can be worked around. 
However, the 6.6 pc size of the images used here falls short of 
the size scales resolved in distant extragalactic systems by a fac- 
tor of; 10. 

With a luminosity in excess of 10 7 Lq 
dWard-Thompson & Robsonl 1 19901) and a mass >10 7 M©, 
W49A is one of the most active star-for ming regions in th e 
Galactic disk, at a distance of 11.4 kpc dGwinn et alJ 1 1992h . 
The region is compose d of a number of optically obscured, 
compact HII regions dWilner et al.l 1200 lb surrounded by a 
molecular cloud w ith a mass >10 5 M© ( Muf son & Liszill 19771 
Simo n et all 1200 ll) . Spectrally, W49A is extremely complex, 
containing numerous features contributed by W49A itself, as 
well as additional clouds associated with t he Sagittarius spira l 
arm which crosses the line-of-sight twice (iPlume et al.l 12004). 
Fig. Q] shows the structure of W49A at infra-red vs sub-mm 
wa velengths and at spat ial scales from ~30 to <1 pc. 

Welch et al. (1987) estimate the cloud density to be 
>3xl0 5 cirT 3 within the inner 1 pc (the core region) and 
~ 10 4 cm 3 in the outer parts of the cloud. Base d on observations 
of the CO J=7-6 transition, iJaffe et al.l (Il987h estimate kinetic 
temperatures of >70 K in the centre, falling below 50 K in the 
outer parts of the cloud. From sub- millimetre continuum maps 
of the dust emission towards W49A. IWard-Thompson & Robsonl 
(1990) find a dust temperature of 50 K and a dust mass 2400 Mq 
over an effective emitting region of ~2 pc. From near-infra- 
red IRAS data they infer an additional population of small, hot 
(T ~350K) dust grains or Polycyclic Aromatic Hydrocarbons 
(PAHs). 

Before W49A can be used as an extragalactic template, a 
better understanding of its physical and chemical properties 
is needed. Rotational transitions of interstellar molecules help 
reaching this goal in two ways: they probe the current density, 
temperature, and optical depth of the gas and they measure the 
abundances of molecular species, which can be compared with 
chemical models to study the evolution of the chemical inventory 
and the physical conditions. Two specific tracers of the chemical 
state of molecular gas are the HCN/HNC ratio and the deutera- 
tion fraction (D/H ratio) of the gas. 

The HCN/HNC abu ndance ratio has long be en known to vary 
with temperature dGoldsmith et alJll98H,[l986l) . Although HNC 
and HCN are assumed to be formed at the same rate (through 
dissociative recombination of the HCNH + ion), the reactions 
of HNC with O and with H, which destroy HNC, have energy 
barriers which make them very slow at low temperatures. As a 
consequence, the HC N/HNC ratio increases with increasing 7kj n 
dSchilke et alj[l992b . 

Observations of molecular D/H ratios (i.e. the abundance of 
a chemical species where one or more D atoms have been sub- 



Table 1. Summary of the observations. 



Line 


Transition 


Frequency 




rms (7 mb ) 






GHz 


S 


K 


HC I8 + 


J=4-3 


340.631 


540 


0.08 


HC I5 N 


J=4-3 


344.200 


120 


0.19 


H 13 CN 


J=4— 3 


345.340 


480 


0.09 


H 13 CO + 


J=4-3 


346.998 


480 


0.12 


HCN 


J=4-3 


354.505 


480 


0.07 


HCO+ 


J=4-3 


356.734 


240 


0.12 


DCO + 


J=5-4 


360.170 


480 


0.07 


DCN 


J=5-4 


362.046 


480 


0.06 


HNC 


J=4-3 


362.630 


583 


0.03 


H 13 CO + 


J=3-2 


260.255 


38 


0.22 


HCN 


J=3-2 


265.886 


19 


0.52 



" : Integration time per pixel 

h : Noise level measured over the central l'xl.5' 



stituted for the H atoms, relative to the abundance of the "nor- 
mal" hydrogen-bearing species) also probe the kinetic temper- 
ature, as well as the chemical his tory, of the ga s. Although the 
cosmic D/H ratio is only 3xl0~ 5 (Linskv 2006), chemical frac- 
tionatio n concentrates deuterium in mol e cules at temperatures 
<50 K dWatsonll 19761 iMillar et al]fl989l) . lHatchell et al] (l998) 
measured DCN/HCN ratios towards a number of hot molecular 
cores (7\i n > 100 K), associated with massive protostars, and 
found ratios of a few times 10~ 3 . These results, and other molec- 
ular D/H ratios in warm gas, can be explained by assuming that 
the deuterium fractionation was set at temperatures <50 K and 
preserved in the icy mantles of interstellar grains, until the re- 
gion was heated and the ices evaporated into the gas phase (se e 
also lLoinard et al Il2002t IParise et al.l 12002b IVastel et al.l l2003). 
Rodg ers & lyjiilarl 01996) showed that deuterium fractionation 
of neutral species can survive for 1 4 — 1 s years in hot gas, 
post-evaporation. Deuterated molecular ions (e.g. DCO + ), on the 
other hand, are destroy ed much more rapidly when the gas is 
warmed to >30 K (e.g. iRoberts & Milldl2007l) and are not ex- 
pected to be abundant in protostellar regions. 

2. Observations 

This work is based o n data obtained as part of the Spectral 
Legacy Survey (SLS; Plu me et al ] 120071) . being conducted at 
the James Clerk Maxwell Telescope (JCMT) on Mauna Kea, 
Hawai'i. The SLS is performing spectral imaging of 2' x2' fields 
towards four targets representing different star-forming environ- 
ments. Once completed, the spectral range covered by the SLS 
will be 330-362 GHz. This paper presents observations of HCN, 
HNC, HCO+, and isotopomers towards W49A (TableG}. For this 
data set, the "central position" refers to R.A. = 19:10:13.4; Dec. 
= 09:06:14 in J2000 co-ordinates. 

The SLS uses the 16- pixel HARP receiv er (325-375 GHz) 
and the ACSIS correlator (buckle et aT1l2009l) . The observations 
were carried out using the HARP4_mc jiggle position switch 
mode, which produces maps sampled every 7.5", and the spectra 
were calibrated by observing at an off-position 14' to the north- 
east of the source. The angular resolution of the JCMT is ~15" at 
345 GHz, which is ~0.8pc at the distance of W49. The spectral 
resolution of th ese observations is ~0.8 km s and the beam ef- 
ficiency is 0.6 dBuckle et al.l F2009). Pointing was checked every 
hour and the pointing is estimated to be accurate to 1.5". 
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Fig. 1. Images of the W49A region at different wavelengths and spatial scales. Left: Three colour MSX image of W49, showing 
the clustered nature of sources in the region: 8.28/vm (red); 12.13/im (green); 14.65/im (blue). Centre: Dust continuum emission 
at 850/mi (SCUBA archive data; see text). Contours are at 20, 40, 60 and 80% of the peak inten sity. This map is ~ 1.5x2' in size, 
slightly smaller than our HARP maps. Right: BIMA observations of hot cores, traced by CH3CN dWilner et al.ll2001l) . We note that 
the central pixels of the HARP maps cover a cluster of several hot cores. 



The data were reduced using the ORAC Data Reduction 
pipeline (ORAC-DR). This automatically processes time-series 
ACSIS data, applying reduction recipes which are determined by 
the incoming data type, and, for HARP, outputting gridded data 
cubes. For SLS data, these recipes include: checking for consis- 
tency between the calculated r sys during the observations and 
the measured r.m.s. noise; checking for variations in the r.m.s. 
noise measured by each receptor across the map; checking for 
r.m.s. uniformity across the spectral range; removing baselines 
from every spectrum; and co-adding cubes taken at the same po- 
sition and frequency. If any receptors are particularly noisy, or 
suffer from bad baselines, they are masked out of the time series 
data and not used in creating the final cubes. 



Supplementary observations of selected molecular transi- 
tions have also been imaged by the SLS team using the JCMT 
Receiver A (21 1-276 GHz; beamwidth 20")- These include the 
HCN and H 13 CO + J=3-2 lines, which are presented below. Half- 
beamwidth spaced raster maps were made, covering approxi- 
mately the same field of view as HARP. The 230 GHz data were 
reduced and calibrated with standard Starlink procedures. Linear 
baselines were removed and the data were Hanning smoothed. 
15 pixel x 15 pixel spectral grids of the data were generated for 
each sight-line such that each pixel corresponds to an (8x8)" 
portion of sky. 



We also made SCUBA maps of W49A, based on data from 
the JCMT archive. These data were reduced using the ORAC- 
DR package which calculates the gain from photometry cali- 
brations and the optical depth from sky-dip measurements. The 
850^/m emission map is shown in Fig. Q] Finally we make use 
of maps of W49A in the HCN, HNC and HCO + 7=1-0 lines, 
made with the IRAM 30m telescope by Peng & Wyrowski (priv. 
comm.). 



Table 2. Observed line ratios 



Species 


Ratio 






1=1-0° 


J=4-3 


HCN/HCO+ 


0.9 


0.56 


HCN/HNC 


1.0 


3.7 


Species 


Ratio J=4-3/J= 1-0 




HCN 


0.73 




HCO + 


1.25 




HNC 


0.22 





Notes. Ratios of velocity-integrated line intensities between species in 
the J= 1-0 and J=4-3 transitions (upper table) and between the J=4—3 
and J= 1—0 transition for each of the main species (lower table). The in- 
tegrated line intensities have been measured in a 67" diameter aperture 
centred on the peak of the emission. 
(n) IRAM 30m data from Peng & Wyrowski, priv. comm. 



3. Results 

3.1. Spatial distributions 

Figures [2] and [3] show maps of the integrated intensities of the 
HCN, HCO + , HNC, H 13 CN, H 13 CO + and HN 13 C J=4-3 and the 
HCN J=3-2 lines around W49A. The same source structure, ex- 
tended in the east-west direction, is seen in both the dust contin- 
uum emission (Fig. [TJ and the HCN and HCO + lines, and also in 
HNC but less pronounced. 

Emission from the HCN J=3-2 line is the most extended, 
being seen across most of the (100x100)" map. The HCN and 
HCO + J=4-3 emission extends over (82.5x60)", based on the 
10% contour, and there is a slightly fainter clump of emission to 
the south-west, seen in both lines. The H 13 CO + J=4-3 emission 
extends over (46x38)", while H 13 CN 4-3 is slightly more com- 
pact. The weaker isotopic lines are seen only towards the central 
pixels of the map: HC 15 N and HC ls O J=4-3 are detected in a 
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30 







30 



30 -30 



Aa (Arcseconds) Aa (Arcseconds) 

Fig. 2. Maps of the integrated intensity ( J T* A dV) across W49: HCN 4-3 (top-left); H 13 CN 4-3 (top-right); HCO + 4-3 (middle-left); 
H 13 CO + 4-3 (middle-right); HNC 4-3 (bottom-left); and HN 13 C 4-3 (bottom-right) line intensities. The main lines are summed 
over 150 km s _1 , while the 13 C-substituted lines are summed over 60 km s _1 . The contours show 5, 10, 15, 20, 40 ... 80% of the 
peak integrated intensity, which is 161 K km s _1 for HCN, 72 K km s _1 for H I3 CN, 319 K km s" 1 for HCO + , and 34 K km s _1 for 
H 13 CO + , 52 for HNC, and 4 for HN 13 C. 



(23x23)" region and DCN J=5-4 over (37x23)". The DCO + 
J=5-4 line was not detected. All sizes refer to the 10% intensity 
level. 

Observations of external galaxies often have limited signal- 
to-noise and therefore consider ratios of velocity-integrated line 



intensities. Figure shows the HCN/HCO+, H 13 CN/H 13 CO + 
and HCN/HNC ratios across the W49A map, based on the in- 
tegrated intensities of the J=4-3 lines shown in Fig. [2] with cut- 
offs applied to only bring out significant structure. Table |2]gives 
the values of these ratios integrated over our maps. Over most 
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Fig. 3. Maps of the integrated intensity ( J T* A dV) across W49: H 13 CO + 3-2 (left) and HCN 3-2 (right). The contours show to 100% 
of the peak integrated intensity (in steps of 15%), which is 24 K km s _1 for H 13 CO + and 172 for HCN. 



of the map, the HCN 4-3 line is weaker than the HCO + line, 
giving integrated intensity ratios <1. The situation is reversed 
for the 13 C-substituted species: H 13 CN is generally stronger than 
H 13 CO + .This reversal indicates a high optical depth of the main 
isotopic HCN and HCO + lines, as also suggested by the lim- 
ited dynamic range of the maps in Fig. [4] and the fact that the 
HCN/H 13 CN, HNC/HN 13 C and HCO + /H 13 CO + intensity ratios 
are 5-10. These ratios are wel below the isotopic abundance ra- 
tio of 60, which indicates line optical depths of ~10, assuming 
that the 13 C-substituted lines are optically thin. 



3.2. Line profiles 

Figure [5] shows the HCN and HCO + line profiles observed 
towards 5x5 pixels (37.5x37.5 arcsec) centred at R.A. = 
19:10:13.4; Dec. = 09:06:14, where the molecular emission 
peaks. Many of the profiles are double peaked and the rest are 
asymmetric. All profiles are best fitted by two Gaussians; the 
combined fit profiles are also shown. Fig. [6] shows histograms 
of the number of spectra in each map as a function of the peak 
velocity of each Gaussian fit. The bin width of 3 km s _1 was cho- 
sen to be large enough to obtain significant signal in each bin, yet 
small enough to be able to locate the peaks of the distribution. 
For all three species, the lower-velocity component of each line 
peaks between and +5 km s while the higher- velocity com- 
ponent peaks between +12 and +15 km s _1 . On the other hand, 
our observations of the rarer isotopologues show single-peaked 
line profiles (Figs[9]and[8]) which peak close to the source veloc- 
ity of +7-8 km s 1 (IWelch et al J 1 19871: Uaffe et alJfl987h . Thus, 
the HCN and HCO + lines are self-absorbed and highly optically 
thick. 

In the central region of W49A, the blue-shifted peaks of 
the main isotopic lines are brighter than the red-sh ifted peaks 
by 50-100%, which is a classic sign of infall (e.g. IZhou et alJ 



Il993h . Infall motion has bee n observed toward similar regions: 
on large scales towa rd W43 dMotte et al. 2003) and on smaller 
scales toward W51 Jzhang & Hoi 1 1 997k In the case of W49A, 
iDickel & Auerl (11994 " suggest large scale free-fall collapse of 
the molecular cloud towards a central HII region based on H CO + 
/= 1-0 and 3-2 line profiles. However. lSerabvn et alJ d 19931) find 
double-peaked line profiles in both CS and C S and argue that 
the starburst in W49A is due to the collision of two clouds. Also, 
our maps show that although "infall-type" line profiles are seen 
in the east-west direction, to the north and south there are regions 
where the red-shifted emission is stronger than the blue, sug- 
gesting outflowing motions in this direction (see also § I3.3.4I >. 
Finally, the clump of emission in the south-west corner is char- 
acterised by single-peaked lines close to the expected source ve- 
locity. We conclude that the gas kinematics in W49A is complex, 
with infall and outflow motions probably occurring simultane- 
ously. Imaging of W49A at high spatial and spectral resolution 
is necessary for a better understanding of its gas kinematics. 

In order to investigate the velocity profiles across the maps, 
Fig. [7] shows integrated intensity maps at the velocities of the 
HCN, HCO + , and HNC J=4-3 fit components as well as at the 
expected source velocity. While the spatial distributions are sig- 
nificantly different at the different velocity ranges, the maps of 
all molecules over the same velocity range are very similar. The 
blue-shifted emission is stronger at the east and west edges of 
the map, while the red-shifted emission is more extended in the 
north-south direction. The maps in Fig. [JJ also suggest the ex- 
istence of unresolved substructure: Possible effects of this sub- 
structure will be discussed in § 13.31 

Based on Fig. [JJ we select three regions of the map away 
from the source centre to study in more detail: the "South- West 
clump" (centred at R.A.=19:10:10.6; Dec=09:05:18) is clearly 
seen at the source velocity and in the red-shifted component, 
but is much fainter in the blue-shifted component; the "Northern 
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H 13 CN/H 13 CO + (4-3) 



HCN/HNC (4-3) 




30 







-30 



30 -30 

Aa (Arcseconds) 



30 



-30 



Fig. 4. Maps of the HCN/HCO+, H 13 CN/H 13 CO + , and HCN/HNC J = 4-3 integrated intensity ( J T* A dV) ratios across W49A 
(greyscale and contours). Contours show 10 to 100% of the peak values (in steps of 15%) which are 1.4 for HCN/HCO + , 4.9 
for H I3 CN/H I3 CO + , and 10.3 for HCN/HNC. A cutoff integrated intensity was applied to each transition before calculating the ratio 
to avoid divide-by-zero errors in low signal regions of the maps. The cutoff values are 1.5 K km s _1 for HCN, 0.5 K km s _1 for 
HNC, 3.0 K km s 1 for HCO+, 0.7 K km s 1 for H 13 CN, and 0.3 K km s 1 for H 13 CO + . 



clump" (R.A. = 19: 10:13.6; Dec.=09:06:48) is a peak in the inte- 
grated intensity seen only in the red-shifted component; while 
the "Eastern tail" (R.A.=19:10:16.6; Dec.=09:05:48) is much 
stronger at the velocity of the blue-shifted component. 



3.3. Column densities 

To derive the column densities of HCN, HNC and HCO + , we 
use the lines of their isotopically substituted species, which have 
a lower optical depth than the main species. Figs.[8]and|9]show 
the observed profiles of these lines, along with single Gaussian 
fits. Although the 13 C-substituted lines are fitted reasonably well, 
these lines are broader than the weaker isotopologues and show 
somewhat asymmetric profiles, so their optical depths may not 
be negligible. This is shown in more detail in Fig. [10] which 
compares the line profiles of the HCO + , HCN, and HNC (and 
isotopologue) spectra towards the central pixel of the W49A 
map. The 13 C-substituted lines do peak at ~8 km s _1 , close 
to the source velocity, but, as well as being asymmetric, the 
H 13 CN J=4-3 line FWHM is -1.4 times that of the HC 15 N J=4- 
3 line. For this reason we have, wherever possible, based our 
calculations of the HCO + and HCN column densities and the 

DCN/HCN ratios on the HC 15 N and HC 18 Q + lines. 

We use the RADEX program (IVan der Tak et al l l2007t) to 
estimate column densities from the observed widths and peak 
intensities of the Gaussian fits to the line emission. This pro- 
gram performs statistical equilibrium calculations including col- 
lisional and radiative processes, and treats optical depth effects 
with an escape probability formalism. The program assumes that 
the optical depth is independent of velocity, which is appropri- 
ate since we are modeling the velocity-integrated line intensi- 
ties. Molecular input data for HCN and H CO + are taken from 
the LAMDA database dSchoier et a l. 2005); we assume that col- 
lisional rates do not change upon isotopic substitution which is 



probably valid given the identical electronic structure. We adopt 
the HCN collision rates for HNC, which assumption is a poten- 
tial source of error since the two species have so mewhat differ- 
ent dip ole moments (3.0 vs 3. 3 D). V ery recently. fSarrasin et al.l 
(2010J) and iDumouchel et alj d2010l) have presented collision 
data for HCN and HNC with He; the impact of these data on 
the present work will be analyzed in a forthcoming paper. 

Since the emission is spatially extended, we assume a beam 
filling factor of unity for all lines. Due to the presence of sub- 
structure within the JCMT beam (e.g.. iDickel et al.|[l999h . this 
assumption provides lower limits on the line optical depths and 
molecular column densities. The measured line widths provide 
an estimate of these limits. For example, the width of the H 13 CN 
4-3 line toward the source centre, which is 1 .4x larger than that 
of the HC 15 N line, suggests a 3-5 x higher optical depth than the 
value of r = 0.2 estimated with RADEX. However, since much 
of our analysis is based on line ratios, the optical depths of indi- 
vidual lines are not a major concern. Section |4]provides further 
discussion on the effect of unresolved substructure. 

For HCN and HCO + , the intensities of two transitions are 
known, which allows an estimate of the ambient physical con- 
ditions (kinetic temperature and H2 volume density). In the case 
of HNC, only one transition was observed, and we have adopted 
the temperature and density found for HCN and HCO + . 



3.3.1. The source centre 

We calculate the column density of the central source adopting 
7 kin = 100K and n(H 2 )= 2xl0 6 cm 3 (Tables 3 and 4). Using 
these physical conditions, T ex (H I3 CO + 7=3-2) ~ T^m, and the 
HCO + column density determinations from H 13 CO + 7=3-2, 
H 13 CO + 7=4-3, and HC ls O + 7=4-3 towards the source centre 
all agree to within a factor of 2, with the column density deter- 
mined from the 7=4-3 and 7=3-2 transitions becoming more 
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Fig. 7. Maps of the integrated intensity ( J T* A dV) across W49A over 5 km s _1 velocity ranges. Left column: -1 to +4 kms -1 (the 
blue-shifted peak); middle column: +6 to +11 kms" 1 (the expected source velocity); right column: +12 to +17 kms -1 (the red- 
shifted peak). Top row: HCN J=4-3; centre row: HCO + J=4-3; bottom row: HNC J=4-3. The contours show 10 to 100% of the 
peak integrated intensity in steps of 15%. The peak intensities are 33, 18, and 16 Kkms -1 , for the blue-shifted, central velocity, and 
red-shifted HCN J=43 emission, respectively; 63, 43, and 33 Kkm s" 1 for HCO + ; and 16, 12,and 5 Kkm s" 1 for HNC. 
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Fig. 5. Observed spectral lines (shaded grey) and Gaussian fit 
profiles (black lines) towards a 37.5"x 37.5" grid centred at R.A. 
= 19:10:13.4; Dec: 09:06:14. The X-axis is velocity relative to 
the line rest frequency (km s _1 ); the Y-axis is T m b (K). HCN 
J=3-2 (top); HCN J=4-3 (middle); HCO + J=4-3 (bottom). The 
dashed vertical line on each plot is at the expected source veloc- 
ity of 7 km s . 



discrepant at lower density. Over the range 7\i n = 75 - 200 K 
and n(H 2 )= 5 x 10 5 - 5 x 10 6 cnr 3 (§ 1), the derived column 
density determinations from the H I3 CN 7=4-3 transitions agree 
very well with those from the HC 15 N 7=4-3 tra nsitions, assum- 
ing is otopic ratios appropriate for the local ISM dWilson & Rood! 
11994 . The calculated HCN and DCN column densities fall by a 
factor of 2 as the kinetic temperature increases from 75 to 200 K 
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Fig. 6. Histograms of the peak velocity for the spectral lines ob- 
served towards the centre of W49A (see Fig. [5} measured with 
double Gaussian fits to the line profiles at each map position. 
The number of lower velocity components is black and the num- 
ber of higher velocity components is white. From top to bottom: 
HCN J=3-2; HCN J=4-3; HCO + J=4-3. 



or by almost an order of magnitude if the density increases from 
5 x 10 5 cnr 3 to 5 x 10 6 cm 4 . Over this range of temperature 
and density, the calculated excitation temperatures for the HCN 
and DCN transitions are similar to each other and increase from 
~ 15 K for 7\i n =75 K and n(H 2 )= 5 x 10 5 cnr 3 to almost 40 K 
for 7^= 200 K and n(H 2 )= 5 x 10 6 cnr 3 . 

Tables [3] and |4] list fit parameters and resulting column den- 
sities, optical depths and excitation temperatures calculated by 
RADEX for the central pixels. As we expect, optical depths 
are generally higher for the transitions of the 13 C-substituted 
molecules than for those containing rarer isotopes. The excita- 
tion temperatures are -40 K for the H 13 CO + and HC lg O + J=4-3 
lines (-86 K for the H 13 CO + J=3-2 lines), -16 K for the H 13 CN 
and HC 15 N J=4-3 lines, and - 18 K for the DCN J=5-4 lines. We 
estimate JV(HCO + ) - 10 14 -10 15 cnr 2 , JV(HCN) - a few times 
10 15 cnr 2 , DCN/HCN = 4xl0~ 4 -3xl0- 3 , and HCN/HCO+ = 
3-10. Based on 3 <x r.m.s. noise levels in the DCO + J=5-4 spec- 
tra, we calculate upper limits on the DCO + /HCO + ratios towards 
the centre of W49A to be 3^x1 0~ 4 for 7 ex = 20-40 K. 

We have also estimated the HNC column density towards 
the centre of W49A, from the J=4-3 transition. For 7^=100 K 
and «(H 2 )=2xl0 6 cnT 3 , we find JV(HNC) = 10 14 cnr 2 so the 
HCN/HNC ratio is -50. The result that JV(HCN)/jV(HCO + ) ra- 
tios are > 1 is contrary that found by naively taking line ratios of 
the main 7=4-3 lines (Fig |4j Table [2]), which shows the impor- 
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Table 3. Estimated HCO + column densities towards the source centre. 



Offset" 


Transition 


Velocity 


T b 


FWHM 


T 


T 


N 


yv(Hco + )* 


(",") 




(km s _1 ) 


(K) 


(km s" 1 ) 


(K) 




(10 12 cm" 2 ) 


(10 14 cnr 2 ) 


15.0,-7.5 


H 13 CO + J=4-3 


5.0 


1.2 


10.0 


37.6 


0.04 


4.5 


3.5 


15.0,0.0 


H 13 CO + J=4-3 


4.0 


2.6 


8.4 


38.3 


0.08 


8.3 


6.4 


15.0,7.5 


H 13 CO + J=4-3 


4.0 


2.2 


11.1 


38.1 


0.07 


9.5 


7.3 


15.0,15.0 


H 13 CO + J=4-3 


8.8 


1.4 


14.2 


37.7 


0.04 


7.3 


5.6 


7.5,-15.0 


H 13 CO + J=4-3 


4.2 


0.6 


10.4 


37.3 


0.02 


2.4 


1.9 


7.5,-7.5 


H 13 CO + J=4-3 


4.0 


3.0 


8.7 


38.5 


0.10 


9.9 


7.7 


7.5,0.0 


HC 18 + J=4-3 


4.2 


0.5 


7.4 


38.6 


0.01 


1.3 


7.2 


7.5,7.5 


HC ls O + J=4-3 


4.3 


0.4 


8.1 


38.6 


0.01 


1.3 


7.3 


7.5,15.0 


H 13 CO + J=4-3 


5.5 


1.9 


10.6 


38.0 


0.06 


7.7 


6.0 


0.0,-15.0 


H 13 CO + J=4-3 


3.8 


0.7 


8.2 


37.3 


0.02 


2.1 


1.7 


0.0,-7.5 


HC 18 + J=4-3 


4.4 


0.4 


7.8 


38.6 


0.01 


1.1 


6.0 


0.0,0.0 


HC 18 + J=4-3 


6.2 


0.5 


9.7 


38.7 


0.02 


2.0 


11. 


0.0,7.5 


HC 18 + J=4-3 


6.9 


0.5 


10.2 


38.6 


0.02 


1.9 


11. 


0.0,15.0 


H 13 CO + J=4-3 


5.5 


1.7 


12.3 


37.9 


0.06 


8.0 


6.2 


-7.5,-15.0 


H 13 CO + J=3-2 


6.4 


1.7 


27.3 


86.2 


0.02 


20. 


16. 


-7.5,-7.5 


H 13 CO + J=4-3 


8.8 


1.2 


12.7 


37.6 


0.04 


5.7 


4.4 




11 \_AJ J — 4— J 


O.J 


Z.o 


1 1 s 


^8 A 


O HQ 

U.U7 


1 1 


Q A 


-7.5,7.5 


HC 18 + J=4-3 


8.9 


0.4 


8.6 


38.6 


0.01 


1.2 


7.0 


-7.5,15.0 


H 13 CO + J=4-3 


6.1 


1.0 


13.3 


37.5 


0.03 


4.9 


3.8 


-15.0,-15.0 


H 13 CO + J=4-3 


9.6 


0.5 


5.6 


37.2 


0.02 


1.0 


0.8 


-15.0,-7.5 


H 13 CO + J=4-3 


9.4 


0.6 


9.6 


37.3 


0.02 


2.2 


1.7 


-15.0,0.0 


H 13 CO + J=4-3 


9.4 


1.0 


11.1 


37.5 


0.03 


4.0 


3.1 


-15.0,7.5 


H 13 CO + J=4-3 


7.4 


1.0 


12.0 


37.5 


0.03 


4.4 


3.4 


-15.0,15.0 


H 13 CO + J=4-3 


6.4 


0.6 


8.4 


37.3 


0.02 


1.8 


1.4 



Notes. HCO + column densities are calculated with RADEX assuming 7\j„= 100 K and n(W^)= 2x 10 6 cirT 3 . For the offset positions where HC 18 + 
was not detected, the data for H 13 CO + J=4-3 were used, except for (-7.5",-15.0"), where H 13 CO + J=4-3 was not detected and H 13 CO + J=3-2 
was used. For the other positions, where all transitions were detected, data for the HC ls O + J=4-3 transition are listed and we note that the column 
density determinations from the H I3 CO + J=4-3 line agree with these to within +50%. 
(fl) Offsets from R. A.= 19: 10: 13.4, Dec -0 9:06: 14 
<*> based on 12 C/ 13 C = 77; ls O/ 18 = 560 dWilson & Ro"od][T994h . 

tance of taking excitation effects into account when estimating 
column densities or abundance ratios. 

3.3.2. The South-West clump 

As Figs. [8] and [9] show, the strength of the isotope-substituted 
transitions falls rapidly away from the source centre. Only the 
13 C-substituted J=4-3 lines are observed towards the SW-clump. 
Fig. QT| shows the detected lines along with Gaussian fit pro- 
files. The lines peak close to the expected source velocity and 
are well-fitted by single Gaussian profiles. Again, we have es- 
timated column densities using RADEX, this time varying the 
physical parameters over the ranges T^ m =20-100 K and n(H2) 
= 5xl0 3 -10 6 cm 3 . We find that, unless «(H 2 ) > 10 5 cirT 3 , the 
HCN and HCO + column densities towards the clump are at least 
10 times higher than towards the source centre, due to the high 
critical densities of these transitions. As this would imply frac- 
tional abundances at least 2 orders of magnitude higher at the 
clump position, it seems more likely that this gas also has a rel- 
atively high density of 10 s — 10 6 cirT 3 . 

We have tried to constrain the kinetic temperature by look- 
ing at the upper limits on the DCN/HCN and DCO + /HCO + ratios 
and at the HCN/HNC ratio. The upper level energies of the DCN 
and DCO + 5-4 transitions are ~50 K, though, so we would not 
expect these lines to be very strong at lower temperatures. For 
7\ in = 20 K, the upper limit on the DCN/HCN ratio from the ob- 
servations is 0.01 whereas steady-state chemical models predict 
DCN/HCN ~ 0.02 at 20 K (see Section|43]), slight evidence that 
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the kinetic temperature in the S W clump is higher than 20 K. The 
HCN and HNC column densities predicted by RADEX give a 
HCN/HNC ratio ~7, which suggests that this region has a lower 
temperature than the source centre. 

Table [5] therefore, lists the line parameters and calcu- 
lated column densities assuming T\^ n = 40 K and «(H2)= 
5xl0 5 c m" 3 . This temper ature is consistent with that deter- 
mined bv lJaffe et alJ dl987h for the outer parts of the cloud (see 
§ UK althoug h the density is ten times higher than estimated by 
Welc h et alj (ll987l) . We note that the HCN and HCO + column 
densities based on the 13 C substituted isotopologues, assuming 
12 C/ 13 C = 77, are 4-6 times higher than those based on ob- 
servations of the main lines. The upper limits on the column 
densities of the deuterated species do not put very strong con- 
straints on the deuterium fractionation: DCN/HCN <0.01 and 
DCO + /HCO + < 0.07. Based on the main isotopic J=4-3 tran- 
sitions, the HCN/HCO + column density ratio towards the SW 
clump is ~1 (~4 based on the 13 C lines), similar to the values 
calculated for the source centre. 

3.3.3. The Eastern tail 

Figure[T2](left column) shows the observed line profiles towards 
the E-tail. Again, the lines are well-fitted by single Gaussian pro- 
files, but the peak velocity of the lines are close to km s , 
blue-shifted from the central source velocity. As discussed for 
the SW-clump above, we consider a range of physical conditions 
and find that column densities based on the HCN J=4-3 com- 
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Table 4. Estimated HCN and DCN column densities towards the source centre. 
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Trail f i nn 
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1 ex 
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/VYHCNV J W(DCN) 
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(II n\ 

K , > 




( Kill S ) 




(Km s ) 






( 1 u cm ; 


(1U cm ) (XlU ) 


15.0,7.5 


HC 15 N J=4-3 


8.8 


0.4 


18.1 


16.3 


0.04 


7.8 


3.5 


4.8 


15.0,15.0 


HC 15 N J=4-3 


10.7 


0.4 


14.7 


16.3 


0.04 


7.3 


3.3 


5.8 


7 s n n 

/ . J,U.U 


ni_ 1>I J — L r — J 


O.J 


n o 

W.V 


IH-.O 


1 f\ A 


n no 


1 f\ 
10. u 


7 
I .z 


1 U. 


7.5,7.5 


HC 15 N J=4-3 


7.1 


0.8 


14.1 


16.4 


0.08 


13.6 


6.1 


8.4 


7.5,15.0 


H 13 CN J=4-3 


7.1 


1.6 


15.6 


16.5 


0.17 


31.8 


2.4 


4.1 


0.0,0.0 


HC 15 N J=4-3 


10.4 


1.1 


11.7 


16.4 


0.11 


15.1 


6.8 


6.2 


0.0 7.5 


HP 15 N T-4-3 

i i v uy j — " «j 


9. 1 


0.9 


12.2 


16.4 


0.09 


13.0 


5.9 


5.4 


1 5 


H 13 PN T-4~^ 
1 1 v^i^i j — t j 


6 5 


1 7 


1 1 7 


1 6 5 

1 U.J 


0.18 


78 7 

ZO. / 


2 2 


3 6 


7 ? on 


1 1 J — H — _J 


V.J 


? 

Z.W 


1J.U 


1 (\ f\ 
1U.U 


U.ZZ 
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i n 

J . u 


^ 1 
j> . i 


7^7^ 


1 1 V^rN J— H — J 


Q 1 
o. 1 


z.w 


1 S 1 


10.0 


ft 7 1 
U.Z1 


^7 S 
J 1 .0 


7 Q 
Z. V 


1 7 


15.0,7.5 


DCN J=5^ 


5.0 


0.1 


14.0 


18.3 


0.01 


4.5 


1.3 




15.0,15.0 


DCN J=5-4 


10.1 


0.2 


14.2 


18.3 


0.02 


9.1 


2.8 




7.5,0.0 


DCN J=5^ 


4.2 


0.2 


7.8 


18.3 


0.02 


3.2 


0.4 




7.5,7.5 


DCN J=5^ 


5.2 


0.3 


9.5 


18.3 


0.03 


8.3 


1.4 




7.5,15.0 


DCN J=5^ 


6.8 


0.3 


14.5 


18.3 


0.03 


10.0 


4.1 




0.0,0.0 


DCN J=5-4 


9.3 


0.2 


14.4 


18.3 


0.02 


7.9 


1.2 




0.0,7.5 


DCN J=5^ 


7.2 


0.3 


13.1 


18.3 


0.03 


10.8 


1.9 




0.0,15.0 


DCN J=5^ 


5.8 


0.2 


11.3 


18.3 


0.02 


6.8 


3.1 




-7.5,0.0 


DCN J=5^ 


11.5 


0.3 


7.5 


18.3 


0.03 


5.2 


1.8 




-7.5,7.5 


DCN J=5^1 


9.4 


0.2 


10.4 


18.3 


0.02 


5.8 


2.0 





Notes. DCN and HCN column densities are calculated with RADEX assuming r kin =100 K and n(H 2 )= 2xl0 6 cm -3 . For the offset positions where 
HC 15 N was not detected, the data for H 13 CN are shown; for the other positions, where both were detected, data for the HC I5 N J=4— 3 transition 
are listed and we note that the column density determinations from the H 13 CN J=4-3 line agree with these to within ±35%. 

{a) Offsets from R.A.=19:10:13.4, Dec.=0 9:06:14 

(b) based on 12 C/ 13 C = 77; 14 N/ 15 N = 450 dWilson & Roodl[T994l) . 

Table 5. Estimated column densities towards the South- West clump. 



Line 


Vlsr 




FWHM 


r.m.s. 


T 

1 ex 


T 


N 




(kms-») 


(K) 


(kms- 1 ) 


(K) 


(K) 




(10 13 cnr 2 ) 


HCN J=3-2 


9.4 


6.1 


9.5 


0.52 


11 


4.9 


33 


HCN J=4-3 


9.2 


4.4 


7.0 


0.06 


11 


5.7 


45 


H 13 CN J=4-3 


8.5 


0.5 


4.2 


0.07 


8.3 


0.2 


2.1 


DCN J=5-4 








0.04 


10 




<0.5 


HNC J=4-3 


8.8 


1.25 


4.3 


0.05 


8.3 


0.7 


6.6 


HCO+ J=4-3 


9.3 


5.0 


7.3 


0.12 


13 


1.4 


6.6 


H 13 CO + J=4-3 


9.0 


0.8 


4.2 


0.08 


12 


0.2 


0.5 


DCO+ J=5-4 








0.07 


10 




<0.5 



Notes. Column densities are calculated with RADEX assuming rkm=40 K and n(H2)= 5xl0 5 cm" 3 towards the South- West clump at 
R.A. =19:10:10.6; Dec. =09:05: 18. We have based the column density ratios on the main lines of HCN and HCO + since the I3 C-substituted spectra 
are rather noisy. 



pared with the J=3-2 transition and the calculated HCN/HNC as measured for H 13 CN and H 13 CO + , the 3<x upper limit on the 
ratio suggest that this gas is also at a lower temperature than the DCN/HCN ratio is ~3xl(T 3 . 
source centre. 



We therefore use the same physical conditions for the E- 
tail as for the SW-clump: 7\ in = 40 K and «(H 2 )= 5x10 s cm 3 . 
Table [6] lists fit parameters and resulting column densities. The 
HCN/HCO + and HCN/HNC ratios are very similar to those 
found for the SW-clump. The upper limit on the DCO + /HCO + 
ratio is more stringent and constrains 7\j n to >25 K. The 
DCN/HCN ratio may be 0.02-0.04, but the DCN J=5^l detec- 
tion is uncertain (see Fig. IT2T > as the line only covers two spec- 
tral channels and the signal-to-noise ratio is <3. Based on the 
r.m.s. noise in the DCN spectrum, using the integrated inten- 
sity for HCN and assuming a line width of 3 km s for DCN 



3.3.4. The Northern clump 

The line profiles of the strongest lines towards the N-clump 
are double peaked, but in contrast to the source centre, the 
red-shifted peaks are stronger than the blue-shifted peaks (see 
Figure [12J right column). These line profiles are not consistent 
with spherical infall motions. Since the isotopic species peak 
closely in velocity to the redshifted emission, the bulk of the gas 
appears to be at this velocity. The blueshifted peak is probably 
due to a bipolar outflow or other unresolved kinematic substruc- 
ture. 
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Table 6. Estimated column densities towards the Eastern tail. 



Line 


Vlsr 


^mb 


FWHM 


r.m.s. 


T 

1 ex 


T 


N 




(km s _1 ) 


(K) 


(kms- 1 ) 


(K) 


(K) 




(10 13 cnr 2 ) 


HCN 1=3-2 


-0.1 


5.05 


9.7 


0.52 


10 


3.85 


25.8 


HCN J=4-3 


-0.3 


3.55 


6.7 


0.04 


9.5 


4.01 


33.0 


H 13 CN J=4-3 


-0.8 


0.20 


3.4 


0.08 


8.5 


0.07 


0.7 


DCN J=5^ n 


-1.0 


0.18 


1.7 


0.07 


12 


0.04 


1.2 


HNC J=4-3 


-0.7 


2.25 


3.4 


0.08 


8.6 


1.7 


9.9 


HCO+ J=4-3 


-0.6 


5.07 


5.5 


0.16 


13 


1.4 


5.0 


H 13 CO + J=4-3 


-0.9 


0.59 


3.2 


0.10 


12 


0.1 


0.3 


DCO+ J=5-4 








0.06 


>10. 




<0.2 



Notes. Column densities are calculated with RADEX assuming 7\j n = 40 K and n(H2)= 5xl0 5 cirT 3 towards the Eastern tail at R.A.= 19:10:16.6; 
Dec. =09:05:48. We show column density ratios based on both the main lines of HCN and HCO + and the 13 C-substituted lines, since there is some 
excess emission at higher velocities in the main HCN and HCO + lines that has not been fitted (see Fig. l 12b . 
{a> Tentative detection based on the r.m.s. noise in the spectrum, iV(DCN)<4xl0 12 cnr 2 for I" cx > 10 K and DCN/HCN <0.002. 



We have fitted single Gaussian profiles to the detected lines 
that are not strongly self-absorbed (HNC, H 13 CN, and H 13 CO + 
J=4-3 and DCN J=5-4; though the signal-to-noise ratio of the 
DCN line is <2) and find that these peak at 15-16 km s _1 
(-19 km s- 1 for DCN J=5-4). The calculated jV(HCN)/JV(HNC) 
ratio, based on the H 13 CN observation, is ~20, which is between 
the values for the source centre (~50) and the two clumps (5-7). 
We therefore adopt an intermediate value for the kinetic temper- 
ature (see also § 14.21 ). Note that our fit to the HNC line profile is 
dominated by the red-shifted peak and thus underestimates the 
HNC column density, but only by ~10%, which does not change 
our result. 

Table [7] lists the fit parameters and resulting column densi- 
ties for 7\m = 75 K and n(H2) = 5xl0 5 cirT 3 . The calculated 
DCN/HCN ratio is an order of magnitude lower than towards 
the E-tail, while the HCN/HCO + ratio is > 1.5 times higher than 
towards the other regions we looked at. 

4. Discussion 

We have presented velocity-resolved maps of the W49A re- 
gion in HCN, HNC, HCO + , and isotopic line emission over 
a 2' (6.6 pc) field at 15" (0.83 pc) resolution. The HCN and 
HCO + emission extends over 60-100", while HNC, H 13 CN and 
H 13 CO + emission is confined to ~40" and rarer isotopes are only 
seen at the source centre. The maps show a molecular 'core', 
elongated East-West, with 'tails' to the North and East, and 
a separate 'clump' to the South-West. While most of the rare 
isotopic species show single-peaked velocity profiles, the main 
HCN and HCO + species show double-peaked, asymmetric pro- 
files suggestive of infall around the peak of the emission. 

Table [8] summarizes our adopted physical conditions and de- 
rived column densities for the various gas components in the 
W49A region. The HCN/HCO + column density ratio is seen to 
vary from ^3 in the SW clump and E tail to «10 in the core and 
the N clump. The HCN/HNC column density ratio is generally 
higher and varies from ^6 in the E tail to «20 in the other gas 
components. The HCN/HCO + column density ratios are quite 
different from the observed line ratios which range from 0.3 to 
1.4, while the isotopic line ratio H 13 CN/H 13 CO + of 0.3-5.3 is 
closer to the column density ratio. 

The HNC line has a similar central velocity and line width 
as the H 13 CO + and H 13 CN lines, which is considerably different 
from the main HCN and HCO + lines. The HNC line thus seems 
to trace the inner core regions more than HCN and HCO + . This 



Table 8. Summary of adopted physical conditions and derived 
column densities for the various gas components in the W49A 
region 



Component 


7\in 


i(H 2 ) 


W(10 14 cm- : 






(K) 


(10 5 cnT 3 ) 


HCN 


HCO+ 


HNC 


Core 


100 


20 


20-70 


2-10 


«1 


SW clump 


40 


5 


10 


4 


0.7 


Etail 


40 


5 


6 


2 


1 


N clump 


75 


5 


10 


0.7 


0.5 



difference may be due to a low HNC abundance in the warm 
central gas, but a possible difference in critical density between 
HCN and HCO + on one hand a nd HNC on the other may also 
play a role ( Sarrasi n et alJl2010r) . 

Since the HCN, HNC and HCO + main isotopic lines are op- 
tically thick, the observed line intensity is not a direct measure 
of column density, but is rather the product of surface area, fill- 
ing factor and temperature. The cloud structure of W49 is not 
homogeneous, as shown by the difference in line width between 
the main and rare isotopic lines, and by the main beam bright- 
ness of the lines which is well below the estimated kinetic tem- 
perature. A possible geometry, commonly invoked to explain 
the large-scale CO emission from Giant Molecular Clouds, is 
an ensemble of small gas clumps in the beam. The inhomoge- 
neous source structure does not directly affect our estimates of 
the beam-averaged column density, which are based on rare iso- 
topic lines. The column densities are indirectly affected through 
our estimate of the volume density, which influenced the fraction 
of HCN in the J -A state. 

We now discuss our derived column densities toward W49A 
by comparing them to predictions from physical and chemical 
models (§ 14. H and 14.2b and to the results of extragalactic obser- 
vations (§ 14. 3b . 

4.1. Comparison to PDR/XDR models 

One interesting co mparison which these observations allow are 
with the models by Meii erink. Spaans. & Israeli (l2007h who pre- 
dict the molecular line emission from irradiated dense gas in 
galaxies using models of PDRs (photon dominat ed regions) and 
XDRs (X-ray dominated regions). Specifically Meii erink et alJ 
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Table 7. Estimated column densities towards the Northern clump. 



Molecule 


Vlsr 




FWHM 


r.m.s. 


T 

1 ex 


T 


N 




(km s~') 


(K) 


(km s _1 ) 


(K) 


(K) 




(10 12 cm" 2 ) 


H 13 CN J=4-3 


15.9 


0.35 


7.4 


0.06 


10 


0.09 


13.8 


DCN J=5-4 D 


18.8 


0.08 


3.5 


0.05 


14 


0.01 


4.3 


HNC J=4-3 


15.7 


1.31 


5.8 


0.05 


10 


0.38 


46.8 


H 13 CO + J=4-3 


15.3 


0.42 


5.4 


0.07 


15 


0.04 


0.9 


DCO+ J=5-4 








0.06 


15 




<0.7 



Notes. Column densities are calculated with RADEX assuming r km = 75 K and rc(H 2 )= 5xl0 5 cm" 5 towards R.A.=19:10:13.6; Dec.=09:06:48. 
(fl) Tentative detection based on the r.m.s. noise in the spectrum, iV(DCN) <3.7xl0 12 cm" 2 for T cx >10 K and DCN/HCN <0.008 
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Fig. 8. As for Fig.|5]but, top to bottom: H 13 CO + J=3-2; H 13 CO + 
J=4-3; HC ls O + J=4-3. 
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Fig. 9. As for Fig. |5] but, top to bottom: H 13 CN J=4-3; HC 15 N 
J=4-3; DCN + J=5-4. 
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Fig. 10. Line profiles towards the central position 
(R.A.=19:10:13.4; Dec.=9:06:14). Note the difference in 
beam size between the 4-3 and 3-2 transitions. 



predict the intensity ratios in the 7=1-0 and 7=4-3 transitions of 
HCN, HCO + and HNC. The measured velocity-integrated values 
of these ratios towards W49A are summarised in Table [2] where 
the values for the 7=1-0 tra nsition come fro m maps made with 
the IRAM 30m telescope bv lPeng et ail d2007l) and Peng et al (in 
prep). 

The observed HCN/HCO+ ratios in both the J=l-0 and 
J=4-3 transitions are inconsistent with all the iMeiierink et al.l 
XDR models, but are consistent with emission from a PDR. The 
J=l-0 line ratio implies a well constrained number density of 
n ~ 4 x 10 4 cm 4 , consistent with the ratio in the J =4-3 tran- 



molecular gas towards W49A 
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Fig. 11. Observed line profiles (shaded grey) and Gaussian fits 
(black lines) for the transitions observed towards the "south-west 
clump" position: R.A.=19:10:10.6; Dec.=09:05:18 



sitions. However, the PDR model has difficulties explaining our 
observed HNC/HCN ratios. Our observed HNC/HCN ratio in the 
J=l-0 line would imply number densities > 3 x 10 5 cirT 3 , while 
the J=4-3 transitions imply n < 2 x 10 4 cm 4 . The HNC/HCN 
line ratios are better explained by the XDR model, where both 
transitions are consistent with densities between 10 4 cm 4 and 
~ 10 5 cm 4 but require X-ray fields with fluxes > 5 erg cirT 2 s _1 
at the lowest densities up to > 100 erg cm 4 s at the highest 
densities. We conclude that no single PDR or XDR model is able 
to reproduce our observations, so that a more complex, multi- 
component model must be invoked. One caveat is the Meijerink 
et al models sh ould be updated with the new collision data for 
HCN and HNC (ISarrasin et al.ll2010h . 

4.2. Comparison to chemical models 

We also compare the observed column density ratios with pre- 
dictions from kinetic models of interstellar chemistry. The mod- 
els consider 340 species linked by > 10,000 reactions with 
rate co-efficients based on the latest re lease of the UMIST 
Database for Astrochemistry, RATEOffl dWoodall et al.ll2007l) . 
The chemical reaction set has been extended to include 
deute rated analogues of all hydrogen-containing species (see 
e.g. iRobe rts & Millar 2000). The models are "pseudo-time- 
dependent", in that the physical conditions remain constant over 
time and they consider only gas-phase chemical reactions (with 
the exception of H2, HD, and D2 formation on grain-surfaces). 
Fig-EJshows steady-state fractional abundances and abundance 
ratios predicted by a set of models run with n(H2) = 10 6 cm 4 
and values for 7\j n between 10 and 200 K. 



1 http://www.udfa.net 
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Fig. 12. Observed line profiles (shaded grey) and Gaussian fits 
(black lines) for the transitions observed towards the "East tail" 
(left) and the "North clump" (right) positions (R.A.=19:10:16.6; 
Dec.=09:05:48 and R.A. = 19:10:13.6; Dec=09:06:48, respec- 
tively) 
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Fig. 13. Steady-state fractional abundances and abundance ratios 
predicted by gas-phase chemical models run at different temper- 
atures; n(H2) = 10 6 cm 4 . 



As discussed in § \T\ the HCN/HNC ratio increases and the 
molecular D/H ratios decrease as T^n increases. The fractional 
abundances of HCN and HNC fall with increasing T^ n , while 
the abundance of HCO + increases. We expect, therefore, that the 
HCN/HCO + ratio will also trace the gas temperature. 

The HCN/HNC ratios we measure in W49A vary between 
50 at the source centre to ~20 towards the N-clump and 6-7 in 
the SW-clump and E-tail. It is these ratios that we have mainly 
used to constrain the temperatures of these regions as 100, 75, 
and 40 K, respectively. In the future, we plan to use the SLS data 
to put more direct constraints on the gas temperature and density 
across the W49A region. 

The HCN/HCO + ratios we observe, however, are between 3 
and 15. The steady-state chemical models predict HCN>HCO + 
only for 7\i n < 25 K. Presumably, though, the source formed 
from a molecular cloud where the gas-phase abundances had 
equilibrated at a much lower temperature. If we use the steady- 
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Fig. 15. Fractional abundance ratios varying over time from a 
gas-phase chemical models assuming n(H2) = 10 6 cm~ 3 and 
7kin= 100 K. The initial abundances come from a steady-state 
model at 10 K. 



Of course, the chemistry in W49A is unlikely to be this sim- 
ple. For each line of sight we may be looking through multiple 
gas-components with different physical conditions (certainly, at 
the source centre the JCMT beam covers several embedded hot- 
core sources; see Fig. []]). The gas-phase chemical models also 
assume that the density and, therefore, the visual extinction are 
high enough that photo-reactions are unim portant. This is sup- 
ported by the IRAC 8 jum image of W49A dPeng et al.ll2010h not 
showing maxima at any of our 'clumps' and 'tails', which would 
be indicative of a PDR. 



10 100 1000 10000 le+05 le+06 

TIME (YR) 

Fig. 14. Fractional abundance ratios varying over time from gas- 
phase chemical models assuming n(H 2 ) = 10 6 cm 3 . The initial 4.3. Comparison to extragalactic regions 



abundances come from a steady-state model at 10 K; the kinetic 
temperatures are 40 K (top) and 100 K (bottom). 



state abundances at 10 K as the initial abundances for the models 
at higher temperatures it takes >10 5 years for the chemistry to 
reach a new steady-state. This is larger than the expected life- 
times of hot-cores in hig h mass star forming regions (3-5 x 
10 4 vr: IWilneret al.ll2001h . 

Fig. [14] shows the evolution of the HCN/HCO + and 
HCN/HNC abundance ratios from gas which has reached steady 
state at 10 K at time=0 yr and is then heated to 40 or 100 K. 
While the HCN/HCO + ratio maintains its "low-temperature" 
value, >1, for >10 4 yr, the HCN/HNC ratio varies much more 
with time, increasing to ~10x its "high temperature" value be- 
fore decreasing to steady-state. It is not really possible from 
these simple models to put an absolute age on the W49A core, 
as the timescales do depend on the density and on the rate of 
warming of the gas, but they could explain why the observed 
HCN/HNC ratios better probe the current gas temperature, while 
the HCN/HCO + ratios are closer to what we expect in colder gas. 

Similarly, for the DCN/HCN ratio, Fig.[l5]shows that it also 
takes >10 4 yr at 100 K for the fractionation set at low tempera- 
tures to begin to fall (see £Q| . In contrast, the DCO + /HCO + ratio 
resets to a "high-temperature" value within 1000 yr. Thus, the 
DCN/HCN ratios of 4xl0~ 4 -10~ 3 observed towards the centre 
of W49A (and the tentative estimate of DCN/HCN = 0.02-0.04 
in the E-tail at 40 K) are not inconsistent with current higher 
kinetic temperatures as, like the HCN/HCO + ratio, the fraction- 
ation may be preserved from an earlier, colder phase of the core's 
evolution. 



Both HCN and HCO + h ave been advocated as tracers of dense 
gas in external galaxies (Ba an et al. I I2008H G racia-Carpio et al.l 
2006), but to date the majority of extragalactic observations have 
focussed on the J=l-0 transitions of these species. However, due 
to their higher critial densities and higher upper energy levels, 
the higher J transitions are likely to be more efficient probes 
of the high density, warm gas associated with star formation. 
In addition, the comparison with lower energy transitions pro- 
vides a probe of the (range of) excitation conditions in the gas 
(Papadopoulos 2007). A survey of nearby active galaxies in the 
J=l-0 to J=3-2 transitions of both HCN and HCO + shows that 
the line intensity from higher J transitions relative to lower tran- 
sitions increases in increasingly starburst-dominated galaxies, as 
does t he ratio of HCO + to HCN emission in a particular J tran- 
sition dKrips et alJl2008l) . The temperature (50 - 100 K), density 
((0.5 -2) x 1 6 cm" 3 ) and column density ( 1 14 - 1 15 cm" 2 ) prop- 
erties of W 49A on the scales observed here are similar to those 
derived by iKrips et all (l2008i) . indicating that W49A may be a 
valuable analogue for understanding the e mission from othe r 
galaxies. Consistent with the observations of lKrips et al.l d2008). 
the strength of the HCO + J=4-3 transition in W49A suggests 
that this line should be a strong and robust tracer of the dense 
gas in extragalactic star forming regions. One advantage of the 
4-3 line over the 1-0 line is that it is much less affected by ab- 
sorption in low-density foreground layers. 

The 7=5-4 lines of HCN, HCO + and HNC h ave been ob- 
serve d for the high-z galaxy APM 08279+5255 dGuelin et al.l 
120071) . In NGC 253, several isotop ologues of HCN HNC 
and HCO + have also been observed (j Hiittemeister et al. 1995; 
IWagg et al.l2005l;lGarcfa-BuriUo et al.l2006HMartm et alj2 006). 
Estimates of ratios derived from these studies are presented in 
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Table 9. Observed line ratios in selected extragalactic objects and comparison to W49A 



Ratio 


APM 08279 


NGC 253 


Centre 


North 


W49A 
East 


Southwest 


HCN/HCO + 4-3° 


0.98 


1.03 


0.5 


0.6 


0.9 


0.8 


H 13 CN/H 13 CO + 2-1 




1.66-3.72* 


«5 C 




«1.4 r 


«1.4 C 


HCN 4-3° / HNC 1-0 


1.57 


0.03-0.14* 


2.0 


1.0 


0.9 


1.4 


H 13 CN/HN 13 C 2-1 




2.29 










HCN 4-3 / H 13 CN 2-1 




0.42-0.56* 










HCO + 4-3 / H 13 CO + 2-1 




0.90-1.51* 










HCO+ 4-3/HC ls O+ 2-1 




2.60 










DNC 2-1 / HNC 1-0 




0.01-0.05* 











7=5-4 for APM 08279 
*: Values for the two observed velocity components 
'": Estimated from the 4—3 ratio (see text) 



Table [9] These studies make these two very active sources par- 
ticularly relevant to compare to the properties derived for each 
component identified in W49A (§§ [33lti334b . 

The comparison between APM 08279+5255 and the com- 
ponents identified in W49A is limited by the number of ratios 
available for the high-z source. In contrast, the data for NGC 253 
are much more complete than those for W49A. Concentrating on 
those ratios where data exist for all sources, we note the follow- 
ing from the data in Table M First, the HCN/HCO+ 4-3 ratios 
in the East and Southwest regions are similar to those in APM 
08279 and NGC 253, whereas the values in the Centre and North 
regions are lower. The minimum value of this ratio measured in 
the entire W49A region is «0.3. Second, the HCN 4-3 / HNC 
1-0 ratios at all positions in W49A are much higher than mea- 
sured toward NGC 253, while the values toward the Centre and 
North regions bracket the ratio measured toward APM 08279. 
Third, the H 13 CN/ H 13 CO + 2-1 ratios, which are based on the 
observed 4-3 ratios, using RADEX at the temperature and den- 
sity of the region to convert the 4-3 ratio to the 2-1 ratio, bracket 
the value measured toward NGC 253. Interestingly, there is lit- 
tle evidence that the fact that APM 08279 is hosting a quasar is 
making any impact on the emission in our observed lines. 

Based on the similarity of the line ratios and the column den- 
sities, we conclude that W49A appears to be a reasonable model 
for extragalactic star-forming regions. 

5. Conclusions 

We have mapped HCN, HNC, HCO + , and related isotopologues 
in the luminous star-forming region W49A. Line ratios of these 
dense gas tracers are often used to probe the chemical state 
and ionisation processes in external galaxies. We find that the 
HCN/HCO + J=4-3 line ratios are <1 across most of the source, 
but that these lines, particularly HCO + , are highly optically 
thick; the H I3 CN/H I3 CO + J=4-3 line ratios are >1 towards the 
source centre and slightly higher just to the north-west; the 
HCN/HCO+ column density ratios are 3-10. The HCN/HNC 
column density ratio varies from «6 in the Eastern tail to *20 
in the other gas components. 

We also find that, although there are embedded hot cores to- 
wards the centre of W49A which are not resolved in these obser- 
vations, column densities towards this region are consistent with 
an average gas temperature of -100 K and density >10 6 cm -3 . 
There appear to be several cooler (7\i n ~ 40 K) but still rela- 
tively dense (n(H2) > 10 5 cm" 3 ) clumps of gas surrounding, and 
possibly infalling onto, the source centre. 



Compa rison of our ob s ervati ons with the PDR and XDR 
models of iMeiierink et al.1 (|2007) shows contradictory results. 
While the HCN and HCO + emission could be consistent with 
a PDR, this is inconsistent with the HNC to HCN ratio which 
points to emisssion from a XDR. We interpret this discrepancy 
as evidence that irradiation plays only a minor part in the chem- 
istry of the molecular gas in W49A. 

We also find that regions of W49A have significant similar- 
ities to observations of the starburst galaxy NGC 253 and the 
quasar host APM 08279. This agreement seems to imply that 
over scales of a few 100 pc, the line emission from external 
galaxies is not dominated by hot core physics but rather by ex- 
tended warm gas-phase chemistry. The W49A region thus may 
be a good template for starburst galaxies, and the current data 
indicate that irradiation by an AGN has a limited effect on the 
chemistry of the surrounding molecular gas. 

The observations discussed here highlight the potential prob- 
lems in the interpretation of extragalactic systems on the basis of 
a few transitions of a limited set of molecular species. The ex- 
tragalatic star forming regions which such observations probe 
are likely at least as complex and structured as W49A. A better 
understanding of the origin of the emission towards W49A and 
the internal physical and kinematic structure of the source will 
require velocity resolved, high angular resolution observations 
with instruments such as the SMA and ALMA interferometers. 
These observations will be able to address the origin of the ve- 
locity structure seen in the JCMT observations as well as deter- 
mine whether different components of the emission are spatially 
segregated due to differences in excitation or chemistry. Together 
with the analysis of the more complete set of molecular species 
probing other physical and chemical environments which will be 
available i n the full JCMT SLS-obs ervations of W49A and the 
Orion Bar (IVan der Wiel et al.ll2009t) . such studies of W49A can 
provide more insight into the relationship between the molecu- 
lar emission, the structure and the star formation in star forming 
regions in other galaxies. 
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